AN EXTENSION OF BOCHNER'S PROBLEM: 
EXCEPTIONAL INVARIANT SUBSPACES 



DAVID GOMEZ-ULLATE, NIKY KAMRAN, AND ROBERT MILSON 

Abstract. A classical result due to Bochner characterizes the 
classical orthogonal polynomial systems as solutions of a second- 
order eigenvalue equation. We extend Bochner's result by dropping 
the assumption that the first element of the orthogonal polynomial 
sequence be a constant. This approach gives rise to new families of 
complete orthogonal polynomial systems that arise as solutions of 
second-order eigenvalue equations with rational coefficients. The 
results are based on a classification of exceptional polynomial sub- 
spaces of codimension one under projective transformations. 
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1. INTRODUCTION AND STATEMENT OF RESULTS 

A classical question in the theory of linear ordinary differential equations, 
which goes back to E. Heine [10] . and which is at the source of many impor- 
tant developments in the study of orthogonal polynomials, is the following: 
given positive integers m and n and polynomials p(x) and q(x) with 

degp = m + 2, degq = m + l, 

find all the polynomials r(x) of degree m such that the ordinary differential 
equation 

(1) p(x)y" + q(x)y' + r(x)y = 0, 

has a polynomial solution of degree n. If there exists a polynomial r(x) 
solving Heine's problem, then it can be shown [18] that for that choice of 
r(x) the polynomial solution y of (pQ) is unique up to multiplication by a non- 
zero real constant. Furthermore, it can also be shown that given polynomials 
p(x) and q(x) as above, a sharp upper bound for the number of polynomials 
r(x) solving Heine's problem is given by 

(2) <7nm ■ = 

\ n 

A well known interpretation of the bound ([2]) is given through an oscilla- 
tion theorem of Stieltjes [17], which says that if the roots of p(x),q(x) are 
real, distinct, and alternating with each other, then there are exactly a nm 
polynomials r(x) such that ([I]) admits a polynomial solution y of degree n. 
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Furthermore, the n zeroes of these solutions are distributed in all possible 
ways in the m + 1 intervals defined by the m + 2 zeroes of p{x). This result 
also admits a physical interpretation in the context of Van Vleck potentials 
in electrostatics, where the roots of y{x) are thought of as charges located 
at the equilibrium configuration of the corresponding Coulomb system. 

The case m = is an important subcase of the Heine-Stieltjes problem. 
The bound g u q = 1 is exact; equation ([T]) is a variant of the hyper geometric 
equation that recovers the classical orthogonal polynomials as solutions of 
indexed by the degree n. In this context, a related classical question, 
posed and solved by Bochner p], specializes the Heine-Stieltjes equation (pQ) 
to an eigenvalue problem. 

Theorem 1.1 (Bochner). Let 

(3) T(y) = p(x)y" + q(x)y' + r(x)y 

be a second- order differential operator such that the eigenvalue problem 

(4) T{P n ) = \ n P n , 

admits a polynomial solution P n (x), where n = degP n , for every degree 
n = 0,1,2,.... Then, necessarily the eigenvalue equation (J3|) is of Heine- 
Stieltjes type with m = 0; i.e., the coefficients ofT are polynomial in x with 
degp = 2, degg = 1, degr = 0. 

If the above theorem is augmented by the assumption that the sequence of 
polynomials {P n (x)} n >o is orthogonal relative to a positive weight function, 
then the answer to Bochner's question is given precisely by the classical 
orthogonal polynomial systems of Hermite, Laguerre and Jacobi, as proved 
by Lesky, [12]. 

Remark 1. Although the literature in the past decades has referred to 
the result above as Bochner's theorem and still the name of Bochner is 
widely associated with the result, in recent years it has become clear that 
the question was already addressed earlier by E. J. Routh, [16] (see page 
509 of Ismail's book pi]). 

If we consider differential operators (|3|) with rational coefficients, say 

(5) p(x) =p(x)/s(x), q(x) = q(x)/s(x), r(x)/s(x), 

where p, q, f, s are polynomials, then the eigenvalue equation (jlj) is, after 
clearing denominators, just a special form of the Heine-Stieltjes equation 
([]]), namely 

(6) p{x)y" + q{x)y' + (f (x) — Xs(x))y = 0. 

It is therefore natural to inquire whether it is possible to define polyno- 
mial sequences as solutions to the Heine-Stieltjes equations with m > ? 
In the present paper, we show that this is possible by weakening the as- 
sumptions of Bochner's theorem. Namely, we demand that the polynomial 
sequence {P n (x)}^ =m begins with a polynomial of degree m, where m > 
is a fixed natural number, rather than with a constant Pq. If we also impose 
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the condition that the polynomial sequence be complete relative to some 
positive-definite measure, then the answer yields new families of orthogonal 
polynomial systems. 

Let us consider the case m = 1. Let b ^ c be constants, and let 

(7) p(x) = k 2 (x - b) 2 + k^x - b) + k 

be a polynomial of degree 2 or less, satisfying k$ = p(b) ^ 0. Set 

(8) a = l/(b-c), 

(9) q(x) = a(x-c)(k 1 (x-b) + 2k ) 

(10) f(x) = -a(ki(x- b) +2k ) 
and define the second-order operator 

(11) T{y) := p (x)y" + My + M- y, 

x — b x — b 

Observe that with T as above, the eigenvalue equation is equivalent to 
an m = 1 Heine-Stieltjes equation: 

(12) (x - b)p{x)y" + q(x)y' + (r(x) - \{x - b))y = 
We are now ready to state our extension of Bochner's result 

Theorem 1.2. Let T be the operator defined in (jlip . Then, the eigenvalue 
equation (|4|) defines a sequence of polynomials {P n (x)}^ =1 where n = degP n 
for every degree n = 1,2,3,.... Conversely, suppose that T is a second- 
order differential operator such that the eigenvalue equation @ is satisfied 
by polynomials P n (x) for degrees n = 1, 2, 3, . . but not for n = 0. Then, 
up to an additive constant, T has the form (jlip subject to the conditions © 
(fTOj) and p(b) ± 0. 

To put our result into perspective requires a point of view that is, in 
some sense, the opposite of the one taken by Heine. Given a collection of 
polynomials y{x) we ask whether there exists a p{x) and a q(x) such that 
this collection arises as the solution set of a Heine-Stieltjes problem ([I]). Let 

(13) V n (x) = (l,X,...,X n ) 

denote the vector space of univariate polynomials of degree less than or equal 
to n. Let M = C V n denote a /c-dimensional polynomial subspace of 
fixed codimension m = n + 1 — k. Let T>2(M) denote the vector space of 
second order linear differential operators with rational coefficients preserving 
M. The assumption of rational coefficients is not a significant restriction. 
Indeed, if dimM > 3, then Proposition 13.11 below, shows there is no loss 
of generality in assuming that T>2(M) consists of operators with rational 
coefficients. We now arrive at the following key definition. 

Definition 1.3. If T>2(M) <2 T>2(V n ), we will call M an exceptional poly- 
nomial subspace. For brevity, we will denote by X m an exceptional subspace 
of codimension m. 



4 



DAVID GOMEZ-ULLATE, NIKY KAMRAN, AND ROBERT MILSON 



We will see below that the concept of an exceptional subspace is the key 
ingredient that allows us to generalize Bochner's result to a broader setting, 
and to thereby define new sequences of polynomials as solutions of a second- 
order equation ([T]). We now construct explicitly an X\ subspace for the 
operator (fill as a preparation for the proof of Theorem 11.21 With a, b, c 
related by 

(14) a(b - c) = 1 
we have 

(15) T(x - c) = 

(16) T((x - b) 2 ) = (2k 2 + oJfei)(s - bf + 2k a(x - c) 

(17) T((x - b) n ) = (n- l)(nk 2 + aki)(x - b) n 

+ (n(n - 2)ki + 2(n - l)ak ) (x - 6) n_1 
+ n(n-3)k (x -b) n ~ 2 , n>2, 

For n = 1, 2, 3, . . ., let £n ,b C V n denote the following codimension 1 poly- 
nomial subspace: 

(18) £%\x) = (a(x - b) - 1, (x - b) 2 , ...,(x- b) n ) 

(19) =(x-c,(x- b) 2 , ...,(x- b) n ) , if a ^ 0. 
The above calculations show that T leaves invariant the infinite flag 

(20) £f c £f c • • • c £% b c • • • , 

It is for this reason, that the eigenvalue equation (j4]) defines a sequence 

of polynomials P\(x), P2(x), By construction, each P n € , while 

equation (fT7|) gives the eigenvalues: 

(21) A n = (n- l)(nk 2 + ak%), n > 1. 

Since T has rational coefficients, it does not preserve V n . Hence, T E 
T^2{£n b ) but T £ T>2(P n ), and therefore £n' b is an X\ subspace. This obser- 
vation is responsible for the forward part of Theorem 11.21 A key element in 
the proof of the converse implication (which we regard as an extension of 
Bochner's theorem) is the following result, which states that there is essen- 
tially one X\ space up to projective equivalence. 

Theorem 1.4. Let M C V n be an X\ subspace. If n > 5, then M is 
projectively equivalent to 

x n ) . 

The answer appears to be much more restrictive than one would have ex- 
pected a-priori. The notion of projective equivalence of polynomial sub- 
spaces under the action of SL(2, R), also an essential element of the proof, 
will be defined at the beginning of Section [2l We complete the proof of 
Theorem 11.21 in Section 
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One of the most important applications of Bochner's theorem relates to 
the classical orthogonal polynomials. In essence, the theorem states that 
these classical families are the only systems of orthogonal polynomials that 
can be defined as solutions of a second-order eigenvalue problem. However, 
new systems of orthogonal polynomials defined by second-order equations 
arise if we drop the assumption that the orthogonal polynomial system be- 
gins with a constant. 

We are going to introduce two special families of orthogonal polynomials 
that arise from flags of the form En ', n = 1,2,... and that occupy a cen- 
tral position in the analysis of the second order-differential operators that 
preserve codimension one subspaces. The detailed analysis of these poly- 
nomial systems will be postponed to a subsequent publication |9j. Here we 
limit ourselves to the key definitions and to the statement of our main result 
concerning the X\ orthogonal polynomials. 

Let a / /3 be real numbers such that a, (5 > — 1 and such that sgna = 
sgn/3. Set 

(22) a = hf3-a), b=^, c = b+l/a. 

I p — a 

Note that, with the above assumptions, |6| > 1. We define the Jacobi-type 
X\ polynomials P n a '^\x), n = 1,2,... to be the sequence of polynomials 
obtained by orthogonalizing the sequence 

(23) x-c,(x- b) 2 , (x - bf, ...,(x- b) n , . . . 
relative to the positive-definite inner product 

(24) (P, Q) ajp := J 1 (1 - ~^\y X)(L W ) dx, 
and by imposing the normalization condition 

(25) 



Having imposed (I25|) we obtain 

2 _ (a + n)(P + n 

4(a + n- l)(p + n- 1 



(26) WP^W 2 [a + nw + n) 



where 

2 a +P+ 1 T(a + n + l)T(p + n + l) 
1 ' n (a + /? + 2n + l)r(n + l)r(a + (3 + n + l) 

is the orthonormalization constant of P n a,/3 \ the classical Jacobi polynomial 
of degree n. 

Proposition 1.5. Set p(x) = x 2 — 1 and let 

(28) T{y) = (x 2 - l)y" + 2a (j^f) {(x - c)y' - y) , 
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be the operator defined by equation Then, the X% Jacobi polynomi- 

als P^^\x), n > 1 form the solution solution set of the Sturm- Liouville 
problem given by (jH) and boundary conditions 

(29) lim (1 - x) a+1 (y(x) - (x - c)y'{x)) = 0, 

(30) lim (1 + xf +1 (y(x) - (x - c)y'(x)) = 0. 

X— > — 1 + 

The corresponding eigenvalues are 

(31) X n = (n-l)(n + a + (3) 

Likewise, for a > 0, we define the Laguerre-type X\ polynomials to be the 
sequence of polynomials L^\x), n = 1,2, .. . obtained by orthogonalizing 
the sequence 

(32) x + a + l,(x + a) 2 , (x + a) 3 , . . . , (x + a) n , . . . 
relative to the positive-definite inner product 

(33) (P,Q) a := 6 X P(x)Q(x)dx, 

Jo (x + ay 

and normalized so that 

(— l) n x n 

(34) Li°'(x) = -, rr + lower order terms. 

v ' n w (n - 1)! 



The orthonormalization constants are given by 

a + n 
a + n — 1 



(35) llW^-^VVx, 



where 

(36) a - r(a + n + 1) 



n — . 

n\ 

are the orthonormalization constants for (x), the classical Laguerre poly- 
nomial of degree n. 

Proposition 1.6. Set p(x) = —x, a = —1, b = —a and let 

(37) T(y) = -xy" + ^ ((x + a + l)y' - y) 

x + a 

be the operator defined by (|11|) . Then, the X\-Laguerre polynomials L„ 
form the solution set of the Sturm- Liouville problem defined by (jl]) and 
boundary conditions 

(38) lim x a+1 e- x (y(x) - (x - c)y'(x)) = 0, 

(39) lim x a+1 e~ x (y(x) - (x - c)y'(x)) = 0. 
The corresponding eigenvalues are 

(40) A n = n - 1. 
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Remark 2. Note that the weight factors ([Ml) and (|33|) differ from the clas- 
sical weights only by multiplication by a rational function. Ouvarov [20] has 
shown how to relate via determinantal formulas the sequence of polynomials 
obtained by Gram-Schmidt orthogonalization of the sequence {1, x, x 2 , ... } 
with respect to two weights that differ by a rational function (see also Sec- 
tion 2.7 in [H]). This does not mean however that Ouvarov's formulas 
apply to the Xi-Jacobi and Xi-Laguerre polynomials defined above, be- 
cause although the weights differ by a rational function, the two sequences 
to which Gram-Schmidt orthogonalization is applied are different, i.e. they 
are {1, x, x 2 , . . . } for the classical polynomials but (|23p and (|32p for the 
Xi-polynomials. 

Indeed, let {P n }^ =0 be the sequence of polynomials obtained by Gram- 
Schmidt orthogonalization from the sequence {l,x,x 2 , . . .} with respect to 
the scalar product (f33l) . Ouvarov's formulas relate the sequence {P n }%Lo 
with the classical Laguerre polynomials. However, the polynomials {P n }^L 
are semi-classical [14J: they do not satisfy a Sturm-Liouville problem, but 
only a second order differential equation whose coefficients depend explicitly 
on the degree of the polynomial eigenfunction. This is the case in general 
for rational modifications of classical weights and orthogonalization of the 
usual sequence, [E]. By way of contrast, the Xi-Laguerre and Xi-Jacobi 
polynomials are eigenfunctions of a Sturm-Liouville problem as established 
by Propositions 11.51 and 11.61 

Once this important precision has been made, we are now ready to state 
the following theorem, which is proved in [9] 

Theorem 1.7. The Sturm-Liouville problems described in Propositions 1 1 . 5l 
and \1.6\ are self-adjoint with a semi-bounded, pure-point spectrum. Their 
respective eigenfunctions are the X\-Jacobi and X\-Laguerre polynomials 
defined above. Conversely, if all the eigenfunctions of a self-adjoint, pure- 
point Sturm-Liouville problem form a polynomial sequence {P n }^ ! L 1 with 
degP n = n, then up to an affine transformation of the independent variable, 
the set of eigenfunctions is X\-Jacobi, X\-Laguerre or a classical orthogonal 
polynomial system. 

Remark 3. In general, a classical orthogonal polynomial system {P n }%Lo 
is no longer complete if the constant Po is removed from the sequence. How- 
ever, in some very special cases the first few polynomials of the sequence 
(although solutions of the eigenvalue equation) do not belong to the corre- 
sponding L 2 space, while the remaining set is complete, [1]. This happens 
for instance for Laguerre polynomials L®(x) when a = — k is a negative in- 
teger: the truncated sequence {L~ k }^ =k forms a complete orthogonal basis 
of £ 2 ([0,oo), x~ k e- x )}. 

The new polynomial systems described in Theorem 11.71 arise by consider- 
ing the m = 1 case of the Heine-Stieltjes problem. As was noted above, 
this allows us to define a spectral problem based on the flag of exceptional 
codimension 1 subspaces shown in (118p . This, in essence, is the "forward" 
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implication contained in Theorem II .71 The reverse implication follows from 
Theorem \1.2\ but requires additional arguments that characterize the X\ 
Jacobi and Laguerre polynomials as the unique X\ families that form com- 
plete orthogonal polynomial systems^. The proof of this result will be given 
in the following paper in this series [9]. 

Let us also point out that some Xj polynomial sequences can be obtained 
from classical orthogonal polynomials by means of state-adding Darboux 
transformations El gJ3 However, this does not explain the very restrictive 
answer that we have obtained for what appears to be a rather significant 
weakening of the hypotheses in Bochner's classification. Let us also mention 
that sequences of constrained, albeit incomplete, orthogonal polynomials be- 
ginning with a first-degree polynomial have been studied in [5] as projections 
of classical orthogonal polynomials. 



2. THE EQUIVALENCE PROBLEM FOR CODIMENSION ONE 

SUBSPACES 

As a preliminary step to the proof of Theorem 11.41 we describe the natu- 
ral projective action of SL(2,R) on V n and on the vector space of second- 
order operators. Our main objective here is to introduce a covariant for the 
SL(2,R) action that will enable us to classify the codimension one subspaces 
of V n up to projective equivalence. 

The irreducible SL(2,R) representation of interest here is the following 
action, Ph>P, on V n : 

(41) P=( lx + 5) n {PoQ, PeV n , 
where 

(42) x = C (£) = ^±i? aS-P 7 = l 

7X + o 

is a fractional linear transformation. The corresponding transformation law 
for second-order operators is therefore given by: 

(43) T(y) = ( 1 x + 5nT(y)oQ, 
where 

5x-p 
-jx + a 



(44) y(x) = (-7X + a) n y 



Here, as part of the definition of an OPS, we assume that the inner product is 
derived from of a non-singular measure. 

2 The polynomials in question do not satisfy deg P n = n, but rather have deg Pi = 
and degP„ = n for n > 2. We will not consider them here. 



AN EXTENSION OF BOCHNER'S PROBLEM 9 

Correspondingly, the components of the operator undergo the following 
transformation: 

(45) p= ( T x + 5) 4 (po(), 

q = ( 7 £ + S) 2 (q o C) - 2(n - 1)7(73 + 5) 3 (p o (), 

r = (r o Q — 727(7^ + 5)(g o C) + n(n — l)7 2 (7x + £) 2 (p o £). 

For convenience, let us set the notation V = V n and G = SL(2,R). 
Let Q n (V) denote the Grassmann manifold of codimension one subspaces 
of V, and let ¥V = Qi{V) denote n-dimensional projective space. We are 
interested in the equivalence and classification problem for the G-action on 
Q n {y)- The action of G is unimodular, and so there exists a G-invariant n+1 
multivector, which we denote by u G A n+1 V. Thus, we have a G-equivariant 
isomorphism <p : A n V — > 1/*, defined by 

</>(«l A • • • A u n )(u)u> = u\ A U2 A . . . A u n A u, u G V. 

Next, we define a non-degenerate bilinear form 7 : V — > y* by means of 
the following relations 

x j j (—!)■? , ifj + fc = n, 



\ j\ ' A;! / j , otherwise. 

Equivalently, we can write 

(46) 7" 1 = Y,{-l) j ( U \x j ® x 

3=0 

Note that 7 is symmetric if n is even, and skew-symmetric if n is odd. 
Proposition 2.1. T/ie above- defined bilinear form is G-invariant. 
Proof. Observe that 

Sym 2 V = {p{x,y) G R[x,y] : deg x {p) < n, deg y (p) < n}, 

and that the diagonal action of G on Sym 2 V is given by 

' ax + (5 ay + (3^ 
jx + 5 ' 7y + 5 

It is not hard to see that p(x,y) = (y — x) n is an invariant. Indeed, 
p(x,y) = (7* + S) n (<yy + 5) n (4±r ~ = (£ ~ 



p(x, y) = (7X + 6) n {>yy + 8) n p 



jy + 5 72 + 5 



Since, 



n 



(47) ( y - x) n = J2(-ir-i( n \xjy n -i 

3=0 

we see that 7 is invariant by comparing (j46]) and (|47l) , □ 
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Since 7 is invariant, it follows that 7 1 ocfi; K n V — ► V is a G-equivariant 
isomorphism. This isomorphism descends to a G-equivariant isomorphism 
$: Gn{V) ^PV. 

Proposition 2.2. Let M G Gn(V) be a codimension one subspace. Then, 
$(M) = {u G F: 7(u,u) = for all v G M}. 

In other words, if v\, . . . ,v n is a basis of M, we can calculate <3?(M) by 
solving the n linear equations 

l(vj,u) = 0, j = l,...,n 

for the unknown ueV. 

There is another natural way to exhibit the isomorphism between Gn(Y) 
and PV. Let M G G n (V) be a codimension one subspace with basis 



Pi( x ) = ^PijX 1 
3=0 

Let us now form the polynomial 



1, 



,n. 



fpio 
P20 



(48) q M (x) = det 



PnO 



Pll 
P21 



Pnl 
n — \ 

-nx 



Plj 

P2j 



Pnj 



Pin \ 
P2n 

Pnn 

(-i)V 



The following proposition shows that, up to scalar multiple, this polynomial 
characterizes M. 

Proposition 2.3. With as above, we have <&(M) = (<?m)- 

Henceforth, we will refer to the subspace of V n spanned by qu as the fun- 
damental covariant of the codimension one subspace M C V n . Thanks to 
the G-equivariant isomorphism between codimension one polynomial sub- 
spaces M and degree n polynomials, we are able to classify the former by 
considering the corresponding equivalence problem for degree n polynomi- 
als. The latter classification problem can be fully solved by means of root 
normalization, as one would expect. 

Recall that a projective transformation (142)) is fully determined by the 
choice of images of 0, 1, 00. Therefore, a polynomial can be put into normal 
form by transforming the root of highest multiplicity to infinity, the root 
of the next highest multiplicity to zero, and the root of the third highest 
multiplicity to 1. 

Proposition 2.4. Every polynomial of degree n is projectively equivalent to 
a polynomial of the form 

k 



(49) 



j=2 
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where rj ^ 0,1, j > 2 and where 

n = rioo + n + ni + n 2 + . ■ ■ n k , > n > n\ > n 2 > ■ ■ ■ 

is an ordered partition of n. The signature partition Woo, no, Hi, . . . and the 
roots rj are invariants that fully solve the equivalence problem. 

Note that in (JIS} there is no factor corresponding to the multiplicity ri^] 
the missing factor corresponds to the root at infinity. 

It is instructive at this stage to show the expression of the covariant (qm) 
of various codimension one subspaces M C V n . 

(1) Consider M\ = {l,x, . . . ,x n ) = V n -\{x). The fundamental co- 
variant is qM±(x) = 1. In this case, qM is equal to its own normal 
form; there is a single root of multiplicity n at oo. 

(2) Consider the exceptional monomial subspace: 

M 2 = (l,x 2 ,x 3 ,...,x n ) =£%°{x). 

Section [3] has more details on this example; see equation (I6ip . The 
covariant in this case is qM 2 { x ) = x n ~ l . The normal form of qM 2 ( x ) 
is a; € V n (x); there is a root of multiplicity n — 1 at oo and a simple 
root at 0. 

(3) Consider the subspace 

M 3 = (l,x,x 2 ,...,x n - 2 ,x n ) = £°(x); 

see equation (|62j) . In this case case qM 3 (x) = x. Therefore, M 2 is 
projectively equivalent to M3. In section [3l below, we show that 
both M 2 and M3 are X\ exceptional subspaces. 

(4) Consider a single gap monomial subspace, 

M 4 = {l,x,...,x j - l ,x j+l ,...,x n ). 

In this case, qM 4 (x) = x n ~ ] . Here, the covariant has one root of 
multiplicity j and another root of multiplicity n — j. 

In the next proposition, we classify the codimension 1 subspaces M C V n 
directly, by exhibiting a normalized basis based on the multiplicity of the 
root at infinity. 

Proposition 2.5. Let M C V n be a codimension one polynomial subspace 
such that qM{x) has a root of multiplicity A at infinity and a root of multi- 
plicity fi at zero; i.e., degojv/ = n — A and [i is the largest integer for which 
x^ divides quix). The following monomials and binomials constitute a basis 
ofM: 

(50) {^ + /3;* A }^ + i {■■■'}'; 

Proof. Observe that qM(x) has a root of multiplicity A at infinity and a root 
of multiplicity \i at zero if and only if, up to a scalar multiple, 

(51) q M {x) = (-lf( n \ X --X _ g { -iy( n \ PjX n-0_ 
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A straightforward calculation then shows that 

i(qm,p) = o, 

where p(x) ranges over the the monomials and binomials in (|50|) . □ 

3. OPERATORS PRESERVING POLYNOMIAL SUBSPACES 

As was noted above, the standard (n + l)-dimensional irreducible rep- 
resentation of SL(2,R) can be realized by means of fractional linear trans- 
formations, as per (|4ip . The corresponding infinitesimal generators of the 
sl(2, R) Lie algebra are given by the following first order operators 

(52) T_ = D x , T = xD x - -, T + = x 2 D x - nx. 

A direct calculation shows that the above operators leave invariant T , n (x), 
and are closed with respect to the Lie bracket: 

(53) [T , T±] = ±T±, [T-,T+] = 2T . 

Since sl(2,R) acts irreducibly on V n , Burnside's Theorem ensures that a 
second order operator T preserves V n if and only if it is a quadratic element 
of the enveloping algebra of the s[(2,R) operators shown in (|52p . Thus, the 
most general second order differential operator T that preserves V n can be 
written as 

(54) T = c a T i T 3+ 2 biTu 

i,j=±,0 i=±,0 

where q,- = c,j, 6j are real constants. For this reason, an operator that 
preserves V n {z) is often referred to as a Lie-algebraic operator. 

For the sake of concreteness we formulate results about invariant poly- 
nomial subspaces by assuming that all operators have rational coefficients. 
However, as the following result will show, this assumption does not entail 
a loss of generality. 

Proposition 3.1. Let T be a second-order differential operator as per ([3]). 
Suppose that Pi(x),Qi(x), i = 1,2,3 are polynomials such that Pi,P%,P3 
are linearly independent and such that T{Pi) = Qi. Then, necessarily the 
coefficients p(x),q(x),r(x) of T are rational functions. 

Proof. By assumption, 

(P<> P> PA fp\ / QA 

\P2 P 2 P2 U = \Q2 ), 



and the matrix on the left is non-singular. Inverting this matrix, we obtain 
rational expressions for p, q, r. □ 
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Proposition 3.2. A second-order operator T preserves V n if and only ifT 
is a linear combination of the following nine operators: 

(55) x 4 D xx -2(n-l)x 3 D x + n(n-l)x 2 , 

(56) x 3 D xx - 2(n - l)x 2 D x + n(n - l)x, 

(57) x D XX i xD xx , Dxxi 

(58) x 2 D x — nx, 

(59) xD x , D x ,l 



A proof can be given based on Burnside's theorem and (J54J). For another 
proof, see Proposition 3.4 of [8]. 

Let us observe that Burnside's Theorem does not apply to general poly- 
nomial subspaces M € V n , and therefore for a general subspace M, there 
is no reason a priori for an operator T £ T>z(M) to also preserve V n . In 
addition to (fl8|) . let us define the following codimension 1 subspaces: 

(60) £»(z) = (l,x,x 2 ,...,x n - 2 ,x n -ax* 1 - 1 ) . 

Indeed, an analysis |13^ [7] of polynomial subspaces spanned by monomials 
brought to light two special subspaces: 

(61) £° n >°(x) = (l,x 2 ,...,x n ), 

(62) £°(x) = (l,x,x 2 ,...,x n - 2 ,x n ). 
These two subspaces are SL(2, M)-equivalent, since 

£° n (x) = x n £^{-l/x). 
Extending the analysis beyond monomials we have the following 

Proposition 3.3. The subspaces £n b ,£n, as defined in (I18p and (|60|) . are 

all projectively equivalent. 

For the proof, see Proposition 4.3 of [Sj. Next, we show that all of the above 
subspaces are Xi, that is exceptional invariant subspaces of codimension 
one. 

Proposition 3.4. A basis of T>2(£n b ) is given by the following seven oper- 

b) 4 D xx - 2(n - l)(x - bfD x + n{n - l){x - b) 2 , 
bfD xx -(n-l)(x-b) 2 D x , 
b) 2 D xx , 

b)D xx + (a(x - b) - 1) D x , 

1 \ „ 2o 
+ 2 ( a r ) D x 



ators: 






(63) 


Ji = 


(x — 


(64) 


J2 = 


(x — 


(65) 


Js = 


[x - 


(66) 


Ja = 


(x- 


(67) 


h = 


D xx 


(68) 




(x - 


(69) 


h = 


1. 



b J x — b 

n) D x — an(x — b), 
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The proof is given in Proposition 4.10 of [8]. 

Observe that J5 is an operator with rational coefficients. Hence, J5 pre- 
serves Sn' b , but does not preserve V n . Therefore, £n b is an X\ subspace. 
Because of projective equivalence, so is Indeed, Theorem 11.41 asserts that 
£n ,b and £ are the only codimension one exceptional subspaces. We prove 
this theorem below. In Section 3, we use Theorem 11.41 to establish Theorem 
11.21 our extension of Bochner's theorem. 

4. PROOF OF THEOREM O 

It will be useful to restate Theorem 11.41 in its contrapositive form. 

Theorem 4.1. Let M C V n , n > 5 be a codimension one subspace. If the 
roots of qm(x) have multiplicity less than or equal to n — 2, then, T>2(M) C 

V 2 (V n ). 

In the preceding section, we showed that if M is projectively equivalent 
to £n°, then quix) has one root of multiplicity n — 1 and another root of 
multiplicity 1. On the other hand, if M is projectively equivalent to V n -ii 
then qM has a single root of multiplicity n. Hence, if the roots of qM{x) 
have multiplicity less than or equal to n — 2, then M is not isomorphic to 
£n° nor to V n -i- Theorem 14.11 asserts that, in this case, T>2{M) C T>2{V n ). 
The rest of the present section will be devoted to the proof of this theorem. 

We begin by writing a second-order differential operator with rational 
coefficients using Laurent series: 

00 

T= £ T k 

k=-N 

where 

(70) T k = x k {a k x 2 D xx + b k xD x + c k ), k>-N 

is a second-order operator of degree k, meaning that Tfcfx- 7 ] is a scalar mul- 
tiple of xi +k for all integers j. Henceforth, for a series L(x) = ^ • Lj x J we 
use the notation 

C j (L) = L j . 

Clearly, if T is a differential operator such that T{M) C M, then necessar- 
ily T{M) C V n . The converse, of course is not true. Nonetheless, it is useful 
to first classify all second order operators that map M into V n , because in 
most instances this larger class of operators turns out to preserve all of V n . 
To complete the proof of the theorem, we consider the more restrictive class 
of operators for which T(M) C M for some limited cases. 

The classification of operators T which map M to V n is the subject of 
the subsequent lemmas. Throughout the discussion, we suppose that T 
is a second-order differential operator and M C V n is a codimension one 
subspace such that T(M) C V n - We also suppose that qM(x) has a root of 
multiplicity A at 00, and a root of multiplicity fj, at 0. By Proposition 12.51 
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this is equivalent to the assumption that x J € M for j = 0, . . . , A — 1, and 
j = n - fi+ 1, . . . ,n. 

Lemma 4.2. IfT k is an operator of fixed degree that annihilates three dis- 
tinct monomials, that is if 

T k [x j } = 

for three distinct j, then necessarily T k = 0. 



Proof. Writing T k as in (|70p and applying it to x J gives 

j(j ~ l)afc + jh + c k = 0. 

Since the above equation holds for 3 distinct j, necessarily a k = b k = c k = 
0. □ 

Lemma 4.3. If X > 2, then T k = for all \ k \> n. 

Proof. By assumption, l,x,x n S M. Hence, if \k\ > n the operator T k 
annihilates these monomials, and hence vanishes. □ 

Lemma 4.4. If qM(x) has only simple roots, then T k = for \ k |> n. 

Proof. By assumption, T k [l] = and T k [x n ] = for all \k\ > n. Hence, 

(71) T k = a k (x k+2 D xx + (1 - n)x k+1 D x ), \k\ > n. 

We are assuming fj, = A = 1, and hence, x J + (3jX G M for j = 2, . . . , n — 1. 
This implies that 

C fc+ i(:zy + /3,-x]) = T k _ j+1 [x : >] + PjTklx] = 

for all k > n and all A; < —2, and hence, by (1711) . 

(72) j(j - n)a k - j+ i + (1 - n)/3,a fe = 0, 
for all such j and /j. In particular, for j = 2, we have 

flfc - 1= 2(2^) /52afc ' 



(73) a fc _j = ( ^ — —fi 2 ) a k 



and more generally, 

n - 1 N j 
2(2 -n)' 

for all j < k + 1 — n if k > n, and all j > if A; < —2. 

Let us argue by contradiction and suppose that a k ^ for some k > n or 
for some A; < — n. By ([72]) and ([73]) . we have 

It follows that by setting 

2(2 -n) P2 ' 
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we have, by ([EE]), that 



Qm(x) = -nx 



= — nx 



= —nx 



= — nx 




= —nx(x — r)' 



This contradicts the assumption that all roots of qM( x ) are simple. □ 

Lemma 4.5. Suppose that T k = for k > n. If A < n — 3, then, Tj- = /or 
A; > 3, and 

(74) T 2 = a 2 (x 4 J D a;a; + 2(1 - n^ZX, + n(n - l)x 2 ) 

(75) Ti = a\x s D xx + b\x 2 D x — n((n — l)a± + b\)x 

Proof. By assumption, x n ,x n ~ 1 + P n -ix x ,x n ~ 2 + (5 n -2X X € M; we do not 
exclude the possibility that (3 n -i = or /3 n _2 = 0. For k > 3, 

C k+n ^{T[x n - 1 + /3 n _ix A ]) = nix"- 1 } + p n ^T k+n ^ x [x x ] = 0, 
C k+n _ 2 (T[x n ~ 2 + P n - 2 x X \) = T k [x n ~ 2 } + P n - 2 T k+n -2-x[x X } = 0. 

By assumption n — 1 — A,n — 2 — A > 1. Hence, for = n, by the above 
equations and by assumption, 

T n [x n - 1 ]=T n [x n - 2 }=0. 

As well, 

T k [x n }=0, k>l. 

Hence, T n annihilates three monomials, and therefore vanishes. We repeat 
this argument inductively to conclude that T k = for all k > 3. For k = 2, 
we have 

T 2 [x n - 1 } = 0, T 2 [x n ] = 0, 

and hence T 2 has the form shown in (|74p . Equation (|75p follows from that 
fact that Ti[x n ] = 0. □ 
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Lemma 4.6. Suppose that T k = for k > n. If X = n — 2, then, T k = for 
k > 4, and 

(76) T 3 = a 3 (z 5 A^ + 2(1 - n)x A D x + n(n - l)x 3 ) 

(77) T 2 = a 2 (x 4 D xx + 2(1 - n^ZX, + n(n - l)x 2 )+ 

+ 2/3 n _ia 3 (x 3 Z) a . -nx 2 ) 

(78) Ti = aix 3 Z? xa; + b 1 x 2 D x - n((n - l)a x + &i)x 

Proof. By assumption, x n , x n_1 + /? n _ix n ~ 2 , x n ~ 3 G M; we do not exclude 
the possibility that P n -l = 0. Hence, for k > 4, 

T fe [x n ] = 0, T k [x n - 1 ] + P n -!T k+1 [x n - 2 ] = 0, T k [x™~ 3 ] = 0. 

Since T n+ i = 0, the above relations imply that T n annihilates three mono- 
mials, and hence vanishes. As before, we repeat this argument inductively 
to prove that T k = for all k > 4. For k = 3, we have 

T 3 [x n - 1 } = 0, T 3 [x n ] = 0, 

and hence (|76p holds. As well, 

T 2 [x n ] = 0, T 2 [x n ~ l ] + p n ^T 3 [x n - 2 ] = 0, 

which proves (I77p . Finally, ?i[x n ] = 0, which proves (|78p . □ 

Lemma 4.7. Suppose that T k = /or A; < — n. If X > 3, then T k = /or 
A; < —3, and 



(79) T_ 2 = a_ 2 £> 

(80) T_i = a-ixD-M. + 6_iZ) : 



i A , n 



Proof. By assumption, l,x,x 2 G M. Hence, for all fc < —3 the operator T k 
annihilates 3 monomials, and hence vanishes. Also note that T_2 annihilates 
l,x and that T_i annihilates 1. Equations (|79p (|80p follow. □ 



Lemma 4.8. Suppose that T k = /or < — n. If X = 2 and fi < 2, then 
the conclusions of Lemma \4- 7\ hold. 



Proof. By assumption, l,x G M, and hence 

(81) T fc [l]=0, T k [x} = 0, k<-2. 

As well, x n_/x + Pn-^x 2 G M, with f3 n -^ ^ 0, and hence, for k < -3, 

C fc+2 (T[x™-^ + I3 n - ll x 2 }) = T k+2 _ n+II [x n -^} + P n ^T k [x 2 } = 0. 

If for some particular k < —3 we have that T k+2 _ n+IM = 0, then T k annihi- 
lates l,x,x 2 . Hence, by induction, T k = for all k < —3. □ 

Lemma 4.9. Suppose that T k = /or < — n. If fi = X = 1, then then the 
conclusions of Lemma \4- 1\ hold. 
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Proof. Since A = 1, we have x n + (5 n -\x G M, where (3 n -\ ^ 0. Hence, 
for k < —3, we have 

(82) C fc+2 (T[x n - 1 + l3 n . lX }) = T k+3 _ n [x n ~ l ] + l3 n ^T k+1 [x} = 0, 
Since \i = 1, we have x 2 + /3 2 x G M, and hence, 

(83) C fe+2 (T[x 2 ]) = T fc [x 2 ] + /3 2 T fc+1 [x] = 0. 

Arguing by induction, suppose that for a given k < —3, it has been shown 
that Tj = for all j < k and that T k [x] = 0. Since /9 n _i 7^ 0, ([82]) implies 
that Tfc + i[a;] = 0. Hence, by ([83]) . T/Jx 2 ] = 0, as well. Since 1 G M, we have 

C fc (T[l]) =T k [l] =0. 

Hence, T k = 0. Our inductive hypothesis is certainly true for k = —n, 
and therefore it is true for all k < —3. Furthermore, T_2[x] = 0. Since 
T_ 2 [l] = 0, as well, {(79j) follows. Relation ([80]) follows from the fact that 
T_i annihilates 1. □ 

Proof of Theorem \4-l\ Let M C V n be a codimension 1 subspace with fun- 
damental covariant qM{x). Let T be a second-order operator such that 
T(M) C M. Necessarily, T(M) C V n , and so we can apply the above lem- 
mas. Let A be the maximum of the multiplicities of the roots of qM(X). We 
perform an SL(2, R) transformation (|42p so as to move the root of qM(x) 
with multiplicity A to 00. Since we have assumed that qM has at least 
two distinct roots, we may simultaneously move one of the other roots to 
zero. Thus, without loss of generality, we suppose that 00 and are roots 
of qM(x) with multiplicities A and [i < A < n — 2, respectively, and that the 
multiplicity of all roots of qM{x) is < A. 

Lemmas 14.31 and 14.41 establish that T k = for \k\ > n. Next, we establish 
that T k = for k > 3 and that Ti,T 2 G T>2{T' n ). Here there are two cases 
to consider 

(1) If A < n — 3, then Lemma 14.51 establishes the above claims. 

(2) Suppose that A = n — 2. Then, 1, x, . . . , x n ~ 3 , x n ~ x + (3 n -ix n ~ 2 , x n 
is a basis for M; we do not exclude the possibility n -i = 0. Since 
T[x n ~ 4 ] G M, we have 

f3 n . 1 C n . 1 {T[x n - A }) = C n _ 2 (T[x™- 4 ]), 

which, by Lemma 14.61 is equivalent to 

12/? n _ia 3 = 12a 2 - 8/3 n _ia 3 

Since T[x n ~ 5 ] G M, we have 

20a 3 = 0. 

Therefore, T k = for k > 4, by Lemma 14.61 The above arguments 
establish that a 3 = 0. Therefore, by equations (j76j) ([77]) (j78j) . T-* = 
and r 2 ,Ti eV 2 (V n ). 
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Next, Lemmas 021 EES 1121 establish that T k = for k < -3, and that 
T_2,T_i G T>2(V n )- Finally To G T>2^P n ) by inspection. Therefore, 

2 

k=-2 

is a sum of operators that preserve V n and therefore preserves V n itself. □ 

5. PROOF OF THEOREM O 

As it was noted in Section [H the forward implication of Theorem 11.21 is 
established by equations (fT5j) (fT7|). Here we prove the converse. Thus we 
suppose that T is a second-order differential operator with rational coeffi- 
cients such that the eigenvalue equation (jH) has polynomial solutions P n {x) 
of degree n for integers n > 1, but not for n = 0. Set 

M n = (P 1 ,P 2 ,...,P n ), n>l. 

By assumption, each M n is a codimension 1 subspace. By Theorem II. 4[ for 
every n > 5, either T preserves V n , or M n = £n°- 

Suppose that T € T>2{V n ) for some n > 5. By Proposition 13.21 T is a 
linear combination of operators (|55p - (|58p . However, since T also preserves 
M n+ \ and M n+ 2 5 and since the operators ([55]) ([56]) ff58|) have an explicit 
dependence on n, our operator T must be of the form 

T(y) = p{x)y" + q{x)y' + ry, 

where degp = 2,degg = 1 and where r is a constant. However, such an 
operator satisfies the eigenvalue equation for n = 0, and hence can be 
excluded by assumption. 

Therefore, M n = £n° for all n > 5. Proposition 13. 3^ asserts that for 
n > 5, there exist constants a n , 6 n such that 7kf n is either £ n n ' n or £% n , as 
per ([18]) ([60|) . We can rule out the latter possibility, because by assumption, 
M n does not contain any constants. Hence, M n = £n n,bn , where, for the 
same reason, a n ^ 0. Hence, there exist constants b n ,c n such that 

M n = (x - c n , (x - b n ) 2 , ... ,(x- b n ) n ), n > 5. 

However, x — C5 and x — c n are both a multiple of Pi(x), and hence c n = C5. 
Also observe that every polynomial p E M n satisfies 

(c n - b n )p'{b n ) +p(b n ) = 0. 

However, since P\^P2,P^ also satisfy 

(c 5 - h)y\b 5 ) + y(b 5 ) = 0, 

we can apply the above constraint to y(x) = (x — b n ) 2 and y(x) = (x — b n ) 3 
to obtain 

2(C 5 " h)(b 5 - b n ) + (65 " b n ) 2 = 0, 
3(C 5 - & 5 )(fr 5 " M 2 + (&5 " &n) 3 = 0. 
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The above imply that b n = 65 also. Henceforth, let us set b = 65 = b n , 
c = C5 = c n , a = l/(c — b). We have established that for every n, 

M n = S% b (x) = (x-c,(x- b) 2 , ...,(x- b) n ). 

Hence, by Proposition 13.41 T is a linear combination of the operators (|63p 
- ([69]) . Again, operators J\, J 2 , Jq have an explicit dependence on n, and 
hence, up to a choice of additive constant, T must be have the form 

T(y) = (k 2 J3 + hJi + k J 5 - akiJ 7 )(y) 

= (k 2 (x - b) 2 + k 1 {x-b) + k Q )y"+ 

+ a(h + 2k /(x - b))((x - c)y' - y). 

By assumption, T(l) is not a constant. Hence, by setting 

p{x) = k 2 (x - b) 2 + k\(x - b) + k , 

we demonstrate that, up to an additive constant, T has the form (jlip sub- 
ject to the condition p(b) 7^ 0. This establishes the reverse implication of 
Theorem O 
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